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ABSTRACT 

We use the cosmic shear data from the Canada-France-Hawaii Telescope Lensing Survey to place 
constraints on f{R) and Generalized Dilaton models of modified gravity. This is highly compli¬ 
mentary to other probes since the constraints mainly come from the non-linear scales; maximal 
deviations with respects to the General-Relativity -i- ACDM scenario occurs at fc ~ l/iMpc“^. At 
these scales, it becomes necessary to account for known degeneracies with baryon feedback and 
massive neutrinos, hence we place constraints jointly on these three physical effects. To achieve 
this, we formulate these modified gravity theories within a common tomographic parameteriza¬ 
tion, we compute their impact on the clustering properties relative to a GR universe, and propa¬ 
gate the observed modifications into the weak lensing quantity. Confronted against the cosmic 
shear data, we reject the f{R) { I/aq I = 10“^, n = 1} model with more than 99.9% confidence 
interval (Cl) when assuming a ACDM dark matter only model. In the presence of baryonic feed¬ 
back processes and massive neutrinos with total mass up to 0.2eV, the model is disfavoured with 
at least 94% Cl in all different combinations studied. Constraints on the {|/rj(, | = 10“"^, n = 2} 
model are weaker, but nevertheless disfavoured with at least 89% CL We identify several specific 
combinations of neutrino mass, baryon feedback and f{R) or Dilaton gravity models that are ex¬ 
cluded by the current cosmic shear data. Notably, universes with three massless neutrinos and no 
baryon feedback are strongly disfavoured in all modified gravity scenarios studied. These results 
indicate that competitive constraints may be achieved with future cosmic shear data. 

Key words: ; Cosmology- Modified Gravity Theories - Methods: analytical, statistical, numeri¬ 
cal 


1 INTRODUCTION 

Explaining the late-time acceleration of the Universe first reported in 
Riess et al. (1998); Perlmutter et al. (1999) represents a major chal¬ 
lenge in modern cosmology, and current interpretations mostly rely 
on the inclusion of dark energy components and/or modifications 
to the theory of General Relativity (GR). One important difficulty 
encountered in solving this puzzle relates to the fact that, by con¬ 
struction, the background dynamics in viable dark energy and modi¬ 
fied gravity models are almost indistinguishable (Bertschinger 2006; 
Song, Hu & Sawicki 2007; Brax et al. 2008). These two frameworks 
only really decouple when considering the evolution of matter den¬ 
sity fluctuations and of perturbations associated with the metric. In 
addition, there are various ways in which a modification of gravity on 
large scales could account for the apparent acceleration (Clifton et al. 
2013; Joyce et al. 2014). Exploiting this, many observational probes 
based on large scale structure formation have been proposed to test 
theories of modified gravity, including galaxy clustering (Pogosian & 
Silvestri 2008; Oyazu, Lima & Hu 2008), integrated Sachs-Wolfe ef¬ 
fect in the cosmic microwave background (CMB) anisotropies and its 


cross-correlation with galaxy density (Song, Peiris, Hu 2007), cluster 
abundance (Jain & Zhang 2008; Lombriser et al. 2010), peculiar ve¬ 
locities (Li et al. 2000; Johnson et al. 2014, 2015), redshift-space dis¬ 
tortions from spectroscopic surveys (Guzzo et al. 2008; Jennings et 
al. 2012; Asaba et al. 2013), 21cm observations (Hall et al. 2013) and 
weak gravitation lensing (Heavens, Kitching & Verde 2007; Schimdt 
2008; Tsujikawa & Tatekawa 2008; Simpson et al. 2013; Wilcox et 
al. 2015). 

In this paper we investigate the extent to which current weak 
lensing surveys can constrain departures from GR. In particular, we 
study the signatures of two specific classes of parametrized modified 
gravity theories, the f{R) and the generalized Dilaton models, on the 
cosmic shear measurement extracted from the CFHTLenS (Erben et 
al. 2013). These models are known to cause an enhancement of struc¬ 
ture formation over scales in the range [0.2 - 20] Mpc/i“^, an effect 
which could be detectable with current lensing surveys. In addition, 
the departure of these models from General Relativity is maximal at 
scales of fc ~ l/i“^Mpc, which are difficult to interpret with other 
clustering data due to the large uncertainty in the galaxy bias. This 
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makes the weak tensing approach special, probing modified gravity 
models at the scale of influence of the ‘fifth’ force. 

In its approach, this paper is an extension of Hamois-Deraps, 
van Waerbeke, Viola & Heymans (2015, hereafter HWVH), where 
the same data was used to place joint constraints on baryon feedback 
models and on the sum of neutrino mass. The general idea can be 
understood as follows: on the one hand, the accuracy achieved by 
modern CMB experiments (Hinshaw et al. 2013; Planck Collabo¬ 
ration 2015) on most ACDM parameters is at the percent level; on 
the other hand, the modified gravity effects we are looking for affect 
the baseline signal by up to 20 percent at small scales. It is therefore 
justified to assume a fixed cosmology and search for possible devi¬ 
ations. Any residual uncertainty in the cosmology can thereafter be 
treated as systematic uncertainty in the analysis. While next gener¬ 
ation weak lensing experiments such as RCSLenS^, DES^, KiDS®, 
Euclid^ and LSST® will have enough statistical power to repeat this 
analysis in a full MCMC pipeline, we demonstrate here that we can 
find interesting results with simpler tools and existing data. 

This paper is organized as follows: In §2 we review the theoret¬ 
ical formulation of structure formation in ,f{R) and Dilaton gravity 
theories; §3 describes the theoretical and numerical modelling of the 
weak lensing signal, and details our cosmic shear measurement from 
the CFHTLenS data. In §4 we present and discuss our results, and 
conclude in §5. The baseline cosmological parameters that are used 
throughout our study correspond to the WMAP9 + BAO -l- SN ACDM 
cosmology: h = 0.6898, Dm = 0.2905, Da = 0.7095, Dx = 0, 
w = —1, Us = 0.831 and Us = 0.969. The reason why we did not 
opt for the Planck cosmology is to minimize the effect of the known 
cosmological tension in our model rejection strategy. Otherwise this 
would involve a full MCMC calculation including all cosmological 
parameters and both data sets as in MacCrann et al. (2014), which is 
not necessary in our approach. In the end however, we do marginalize 
over this cosmological discrepancy. 


2 MODIFIED GRAVITY THEORIES 

Modified theories of gravity can be distinguished by their screen¬ 
ing properties in dense environments. Indeed, given the strong So¬ 
lar System constraints, these theories need to have a built-in screen¬ 
ing mechanism, suppressing the deviations from GR. Three types of 
such mechanisms have emerged in the last few years: the Chameleon, 
K-mouflage and Vainshtein models (see Brax & Valageas 2014, 
for a comparison between these different screening mechanisms). 
On the one hand, K-mouflage and Vainshtein models involve non¬ 
linear kinetic terms describing additional scalar fields whose pres¬ 
ence modifies GR predictions. On the other hand, modifications of 
the Chameleon type can be broadly categorized as either contain¬ 
ing additional couplings between the metric and new scalar fields, or 
involving extra geometric terms. These two equivalent descriptions 
can be captured by the tomographic parameterization, which will be 
used throughout this paper (Brax, Davis & Li 2012; Brax, Davis, Li 
& Winther 2012). 

In all Chameleon cases, modifications of gravity induce a global 
enhancement of the effective force of gravity, due to the ‘fifth force’, 
which directly translates into an increase of structure formation. 
In this Section, we review two different types of modified gravity. 


^ http://www.rcslens.org 
^ http://www.darkenergysurvey.org 
® http://kids.strw.leidenuniv.nl 
http://sci.esa.int/euclid 
® http://www.lsst.org/lsst 


namely the Dilaton and the f{R) models; we describe their distinct 
screening mechanisms, and detail their parameterization in the con¬ 
text of large scale structure formation. 


2.1 Gravity in Dilaton models 

The Dilaton and Symmetron^ theories of modified gravity are 
Chameleon models that exhibit the Damour-Polyakov property 
(Damour & Polyakov 1994), according to which the coupling be¬ 
tween the scalar field p and the rest of the matter components 
approaches zero in dense environments (Pietroni 2005; Olive & 
Pospelov 2008; Hinterbichler & Khoury 2010). In contrast to the case 
of f{R) theories, described in Sec. 2.2 below, the scalar field here 
takes on a small mass everywhere and thus mediates a long-range 
(screened) force. 

These Dilaton models are scalar-tensor theories, where the ac¬ 
tion defining the system takes the general form 


S = 


= / 


d^XsJ—g 






( 1 ) 


where Mpi = is the reduced Planck mass (in natural 

units), Aq is the cosmological constant term today, g is the determi¬ 
nant of the Einstein-frame metric tensor and g the determinant of 
the Jordan-frame metric tensor The two metrics are connected 
via a conformal rescaling 

~ 7l {p)g^i/- (2) 


The various matter fields ipm are governed by the Jordan-frame 
Lagrangian density Crn and the scalar field ip by the Einstein-frame 
Lagrangian density tCp = —l/2(Vv3)^ — V{}p), where V{p) is the 
potential of the scalar field^. There is no explicit coupling between 
matter and the scalar fields, and the fifth force on matter particles due 
to p arises from the conformal transformation given by equation (2) 
(more precisely, through gradients of A). 

In the original Dilaton model, the potential V {p) and the cou¬ 
pling® A{ip) with the metric have the following functional forms: 


V(p) = 14 exp 



(3) 


where {V*, ^2 } are the two free parameters. In regions where 0, 
the coupling to matter is negligible, and gravity converges to GR. 
However, the field nevertheless mediates a long range gravitational 
force that has an effect elsewhere, i.e. in less dense environments. 
This model can be generalized to a greater class of Dilaton models, 
by keeping the coupling function as in equation (4) but considering 
more general potentials. Then, instead of specifying the model by 


° We do not further investigate the Symmetron, K-mouflage nor Vainshtein 
models in this paper. 

^ In equation (1), we explicitly added the cosmological constant term Aq, so 
that the minimum of V((p) is zero and is reached for p oo. Alternatively, 
this term could also be interpreted as the non-zero minimum of the scalar field 
potential. 

® This coupling is often defined as A{p) = 1 + ~ where 

-/v2pi 

is some free parameter of the model. We opted to absorb p* into p in 
equation (4), a choice that has no physical impact anyway. 
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its potential it is re-casted in the tomographic parametriza- 

tion {j3[a),m{a)} in terms of the cosmological scale factor a{t), 
where the coupling /3(a) and the scalar field mass m(a) are defined 
as (Brax, Davis & Li 2012; Brax & Valageas 2013): 

^(a) = ^[(^(a)] = (5) 


m {a) = m [ip{a),p{a)\ = 




(6) 


Hereafter, we denote with an overbar unperturbed cosmological 
background quantities, and with a subscript ‘0’ quantities evaluated 
today. For instance, p{a) — Mpi/a^ is the background 

matter density, (p is the mean value of the field, Hq is the current 
value of the Hubbles parameter, and HmO is the current matter den¬ 
sity. Also, c is the speed of light in vacuum. In this paper we consider 
the simple forms 


m(a) = mo a /3(a) 


/3o exp 


’ - 1 


3- 2r 


(7) 


with 

_ 9A2Hmo77o ,0, 

5 9 2 ■ I®) 

c'^mg 

In this framework, the Yukawa potential given by equation (3) corre¬ 
sponds to r = 3/2. The values of the free parameters {mo, r, fio, s} 
that enter equation (7) are displayed in Table 1. The models (A, B, 
C, D} were chosen such as to correspond to those studied in Brax 

6 Valageas (2013) and Brax et al. (2012), where detailed compar¬ 
isons between numerical and analytical calculations are presented. 
More specifically, the models {A, B, C} probe the dependence on 
(s, /3o, n} respectively, all other parameters being fixed, while mod¬ 
els D probe the dependence on mo at fixed A 2 . We added the models 
E that probe the dependence on the parameter mo at fixed (s, /3o, n}. 
These models probe deviations from the ACDM cosmology of less 
than 20%, in terms of the matter power spectrum. Let us now explore 
the detailed mechanism through which the power spectrum of matter 
fluctuations is affected by this theory. 

In these Dilaton models, the coupling function A is always very 
close to unity, so that most Einstein-frame and Jordan-frame quan¬ 
tities (e.g., Hubble expansion rates or densities) are almost identi¬ 
cal. Indeed, using |A — 1| <C 1, we can see from equations (4) 
and (5) that A ~ 1 -|- /3^/(2A2). From equation (8) we also obtain 
A 2 ~ (cmo/Ho)^. Solar System tests of gravity such as that anal¬ 
ysed in Chiba (2003) imply that mo > lO^Ho/c, whence A 2 > 10® 
and 


|A-1|<10“®. (9) 

Therefore, the Jordan-frame and Einstein-frame scale factors and 
background matter densities, related by a = Aa and p = A~*p, 
can be considered equal, as well as the cosmic times and Hubble 
expansion rates. In the rest of this Section we work in the Einstein 
frame, where the analysis of the gravitational dynamics are simpler. 
In the Einstein frame, the Eriedmann equation takes the usual 

form 

SMpiH = p-\-ptp-\-pK, ( 10 ) 

where we explicitly separate contributions from the matter (p) and 
scalar field {p^) components and from the cosmological constant pA. 
The background value of the scalar field potential is given by 


Combining with equation (6), and writing m = 

dtp 3 Pp dV 3P^p^ 

da c^Mpxam? ’ da c^Mp,am? ’ 


whence 



m{a), this leads to 

( 12 ) 

10”®. (13) 


Thus, the scalar field energy density is dominated by its potential 
term, which is negligible as compared with the matter density. There¬ 
fore, the Eriedmann equation (10) is governed by the matter density 
and the cosmological constant and we recover the ACDM cosmolog¬ 
ical expansion, SMpiH^ ~ P + pA, up to an accuracy of 10“®. 

We now briefly consider the behaviour of metric and density 
fluctuations. In the quasi-static limit, the scalar field is given by the 
Klein-Gordon equation. 


c" 2 dV ^ dA 
—V ip = 

dip dp 


(14) 


and at linear order over the matter density and scalar field fluctuations 
we obtain 


5p 

Mp\ 


/3 


5p 


(15) 


c^Afpj m? + tA ja'^ ’ 
where k is the comoving wave number. Using equation (11) this gives 


|M| 


M 

p 


cm 


1 


< 10” 


(16) 


1 -I- 

so that the perturbations of the conformal factor M are negligible 
compared to unity. Also, 

2 1 


Sp 


cm 


I -I- k'^/o?m?- 


< 10“ 


(IV) 


hence fluctuations of the scalar field energy density are negligible 
compared with the matter density fluctuations. 

Therefore, the main source that drives modifications to structure 
growth is not a different background evolution, nor perturbations in 
the scalar field energy density, but really the action of the fifth force 
on the matter field. In the Newtonian gauge, the perturbed metric can 
be written as 

^ 2^)dd + a^(f)(l - 2-i!)&ijdx^dxT (18) 


where 4? and are the Einstein-frame metric gravitational potentials. 
Using equations (15) and (17), we can check that the impact of the 
scalar field fluctuations on the metric potentials are negligible, and 
we have within a 10“® accuracy 


$ = '!' = 


(19) 


where is the Newtonian potential given by the Poisson equation, 

= ^-kQ^Sp = ( 20 ) 

In the above expression, S = 5p/p is the matter density contrast. 
However, the dynamics of matter particles are modified by the scalar 
field, which gives rise to the fifth force given by F = — c^Vln A. 
That is, in the Euler equation we must add a fifth-force potential, 
T'a = c^lnA, that is not negligible. When solving for structure 
growth given the parameters listed in Table I, the new term can lead 
to lO — 20% deviations in the matter density power spectrum. 


2.2 Gravity in f{R) theories 


dp 




In models based on f{R) gravity, the Einstein-Hilbert action is mod¬ 
ified by promoting the Ricci scalar Rto a. function of R (Buchdahl 
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1970; Starobinsky 1980, 2007; Hu & Sawicki 2007). The new action 
S in f{R) gravity theories can be written as: 


^ ~ J d'^x^—g 


Mi 


■[R + f{R)]-At + Cm{^^^) 


( 21 ) 


where we explicitly added the cosmological constant contribution®. 
The f{R) models are most easily described in the Jordan frame, 
which is why, in this Section, we denote with a tilde Einstein-frame 
quantities instead of Jordan-frame ones, contrary to the notation of 
§2.1. In the parameterization of Hu & Sawicki (2007), the functional 
form f{R) can be expressed in the high curvature limit as: 


f{R) 


fRo r ^ f Rq^" 

n i?" ’ dR Rn.+i- 


( 22 ) 


The two independent parameters, and n, can be constrained 
by observations. In the above expression, Rq is the present value 
of the Ricci scalar for the cosmological background. Note that this 
parametrization and that of Starobinsky (2007) both reproduce the 
same results in the large curvature regime. 

The f[R) theories of gravity also invoke the Chameleon mech¬ 
anism to screen modifications of GR in dense environments such as 
in our Solar System. Specifically, this occurs by requiring that all ex¬ 
tra terms vanish in high curvature environment, such that f{\R\ 

I J?o I) —>■ 0. The background expansion otherwise follows the ACDM 
dynamics and the growth of structure is only affected on intermediate 
and quasi-linear scales. 

There is an essential connection between the formulation of the 
f{R) theory presented above and scalar-tensor theories of modified 
gravity. Upon the coordinate rescaling (recall 

that in this Section is the Einstein-frame metric) with A{ifi) = 
exp[/3:/3/Mpi] and f} = Ij x/fi, the f[R) modifications to GR are re- 
casted as arising from contributions of an extra scalar field ip, subject 
to a potential V {p) given by: 


V{p) 


Mgi ( RfR-f{R) \ 

2 V (l + /«)^ J’ 


fn = exp 


2/3y> 

Mpi 


-1. (23) 


In that sense, /(/?) theories are equivalent to a scalar-tensor the¬ 
ory expressed in the Einstein frame (Chiba 2003; Nunez & Soignaik 
2004). In this new formulation, the screening mechanism takes an¬ 
other form: the mass of the scalar field grows with matter density, 
and a Yukawa-like potential suppresses the fifth force in dense en¬ 
vironments. This can be conveniently reformulated by saying that 
screening takes place wherever the scalar field is small compared to 
the ambient Newtonian potential. 

It turns out that all Chameleon-like models such as f{R) theo¬ 
ries can again be parameterized by the value of the mass m(a) and 
the coupling P{a) of the scalar field, in terms of the scale factor a 
and the associated background matter density p(a). With the specific 
functional form of f{R) given by equation (22), we can directly re¬ 
late {n, /flol to {/3(a), m(a)} via: 


m{a) — mo 


( 4Hao + Hmoa ® 
\ 4HaO + 


( n -|- 2)/2 


mo 


Ho / Hmo + 4Hao 


c Y (n+ l)|/flo| 


P{o 


1 

VE' 


(24) 


In this paper, we consider values of n = {1,2} and |/flo| = 
{10“'^, 10“®, 10“®}. The larger value of is currently ruled 


® The terms R and Aq are often included within the function f{R). Written 
in the form of equation (21), f{R) describes deviations from both GR and the 
ACDM cosmology. 


out by other independent probes, so this serves as a consistency test. 
The numerical values for {/3{a), m(a)} corresponding to these three 
models are listed in Table 1. 

As for the Dilaton models described in §2.1, the f{R) models 
that we consider in this paper follow very closely the ACDM cos¬ 
mology at the background level, mainly because | <C 1. Indeed, 
from the action (equation 21) one obtains the Friedmann equation as 
(Tsujikawa 2007); 

3M|i [jT® - + H) + f/6 + IrrMr] =p + pa, (25) 

where the dot denotes the derivative with respect to cosmic time t 
and /rr = f /dR?. In the background we have R = 12H^ + QH 

and we can check that all extra terms in the brackets in equation 
(25) are of order |7T^, so that we recover the ACDM expansion, 
SMpiH^ = p + Pa, up to an accuracy of 10“^ for I/rq \ < 10“'^. 
Moreover, the conformal factor A{p) is given by A = 
so that I A —1| < 10““^ and the background quantities associated with 
the Einstein and Jordan frames can be considered equal (and equal to 
the ACDM reference) up to an accuracy of 10“^. 

Considering the metric and density perturbations, we can again 
write the Newtonian gauge metric as in equation (18) (but this is now 
the Jordan-frame metric). Then, in the small-scale sub-horizon limit 
k/a ';$> H/c, the modified Einstein equations lead to (Tsujikawa & 
Tatekawa 2008) 

+ 47reN(5p, (26) 

V2 

-^4' = -l^^fR + 47rC/N(5p, (27) 

where S/r = Jr — fp and Sp = p — p. Therefore, in terms of the 
Newtonian gravitational potential T'n defined as in GR by equation 
(20), we have 

2 2 

$ = 4/N-y5/fl, V& = 4'N + y5/fl. (28) 

Thus, because we work in the Jordan frame (in contrast with the Dila¬ 
ton case presented in §2.1), the modification of gravity directly ap¬ 
pears through the metric potentials. The fluctuations of the new de¬ 
gree of freedom S/r are given by: 

-2— S/r = SR - SuGnSp. (29) 

Finally, the dynamics of the matter particles is given by the geodesic 
equation, where the Newtonian potential that appears in GR is re¬ 
placed by the potential 4> given in equation (28). 


3 WEAK TENSING 

3.1 Theory 

3.1.1 Weak lensing convergence power spectrum 

In all the cosmologies considered in this paper, we work in the New¬ 
tonian gauge with the perturbed metric given by equation (18), where 
<I> and 4r are the metric gravitational potentials^®. In practice, we 
measure the statistical properties of weak lensing distortions by sum¬ 
ming over many galaxy images. This means that the measured signal 
is an integral over selected sources with a broad redshift distribution 

In the Dilaton models, this is understood as the Einstein-frame metric 
while in the f{R) models this is the Jordan-frame metric, following the ap¬ 
proach described in §2. In any case, we can work in either frame as the obser¬ 
vational results do not depend on this computational choice. 
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Table 1. Parameters describing the modified gravity theories considered in our study, mapped on the {/3(a.), m(a)} surface, parameterized with {mo, r, /3o, s} 
following equation (7). The first five rows correspond to different realizations of the generalized Dilaton theories. The last two rows show f{R) theories with 
n = 1 and 2 respectively, in which mo is given by equation (24), while r = 3{n + 2)/2, 0o = l/-s/6 and s = 0. 


Model 

mo(/i/Mpc) 

r 

00 

s 

(A1,A2, A3) 

(0.334,0.334,0.334) 

(1.0,1.0,1.0) 

(0.5, 0.5, 0.5) 

(0.6,0.24,0.12) 

(B1,B3, B4) 

(0.334,0.334,0.334) 

(1.0,1.0,1.0) 

(0.25,0.75,1.0) 

(0.24,0.24,0.24) 

(C1,C3,C4) 

(0.334,0.334,0.334) 

(1.33,0.67,0.4) 

(0.5, 0.5,0.5) 

(0.24,0.24,0.24) 

(D1,D3, D4) 

(0.667,0.167,0.111) 

(1.0,1.0,1.0) 

(0.5, 0.5,0.5) 

(0.06,0.96,2.16) 

(E1,E3,E4) 

(0.667,0.167,0.111) 

(1.0,1.0,1.0) 

(0.5, 0.5,0.5) 

(0.24,0.24,0.24) 

iz = 1. logio l/flo 1 = (-4, -5, -6) 

(0.042,0.132,0.417) 

(4.5, 4.5,4.5) 

(0.408,0.408,0.408) 

(0,0,0) 

n = 2, logic l/fio 1 = (-4- -5. -6) 

(0.034,0.108,0.340) 

(6.0, 6.0, 6.0) 

(0.408,0.408,0.408) 

(0,0,0) 


n{zs)dzs (mapped to n{x)dx iri terms of the radial distance, given 
the Jacobian dxjdz) that we normalize to unity. Thus, introducing 
the kernel g(x) that defines the radial depth of the survey: 

g{x) = [ dxMXs)— -(30) 

J X 

the integrated convergence field at a position 6 on the sky reads as: 


k{0)=( dx\g{x)'^^^^\{x,xO)- (31) 

Jo 

We assumed a flat background universe in the above equation, and 
introduced the weak lensing potential, defined by 


= 


$ + tp 
2 


(32) 


which is a convenient when computing weak lensing modifications to 
GR. Solving equation (31) in multipole space and taking the ensem¬ 
ble average of the squared complex norm, we obtain the convergence 
power spectrum: 


C7=\ dx^^—P^^dGx\z) (33) 

Jo c X 

as an integral over the weak lensing power spectrum 
Note that the above also assumes both Limber and Born approxi¬ 
mations. From this, we also derive predictions for the cosmic shear 
two-point correlation functions ^±(0), computed as: 

= ^ / CUoifim £ d£ (34) 

where Jo/i{x) are Bessel functions of the first kind. 


3.1.2 Ct in General Relativity 

In the ACDM cosmology -l- GR case, we can exactly express the 
weak lensing convergence power spectrum (33) in terms of the total 
matter power spectrum P{k) via Poisson equation. Indeed, we can 
safely neglect the anisotropic stress, and General Relativity gives: 

$wi = $ = = ^N, (35) 

where is the Newtonian potential given by Poisson equation 
(equation 20). Therefore, we recover 

(36) 

and the convergence power spectrum (33) becomes: 

noo 

C^= dxW{xfP{e/x;z), (37) 

Jo 


with 

W{x)= ^^^°^° g{x)0 + z). (38) 

3.1.3 Ct in theories of modified gravity 

For the Dilaton models, we have seen in equation (19) and in §2.1 
that the two Einstein-frame metric potentials are equal to the Newto¬ 
nian potential up to order 10~® accuracy, and that background cos¬ 
mological quantities such as the Hubble expansion rate and the radial 
comoving distances are equal to those of the ACDM reference within 
that same accuracy. This means that equations (35)-(36) apply as in 
GR, and that Cf is again given by equations (37)-(38). Therefore, in 
terms of this weak lensing statistics, the modification of gravity and 
the departures from the ACDM-l-GR results only appear through the 
modified matter density power spectra P(fc; z), which we describe in 
§3.2.3. 

In the case of /(i?) models, we have seen in equation (28) that 
the two Jordan-frame potentials are different from the Newtonian po¬ 
tentials, receiving contributions from terms linear in 5/ij. However, 
these two extra terms exactly cancel in the weak lensing potential 
(equation 32) such that <E>wi = ^n. Therefore, we recover equa¬ 
tions (36)-(38) in /(i?) models too. Moreover, we have seen that both 
the Jordan-frame and Einstein-frame background quantities are equal 
to the reference ACDM background quantities up to an accuracy of 
10“^ for l/ijo I < 10~^. This means that weak lensing statistics can 
again be computed in the reference background cosmology, so long 
as the modified matter density power spectrum is used. 


3.2 Non-linear matter power spectrum 

The choice of non-linear power spectrum to insert in equation (37) 
depends on the cosmology under investigation. In this paper, we are 
interested in constraining modified gravity models, but with respect 
to a ACDM baseline, these are strongly degenerate with universes 
that include baryon feedbacks and/or massive neutrinos. In the con¬ 
text of cosmic shear, these phenomenas are therefore intrinsically 
connected and must be jointly analysed. We detail in this Section 
how we combine all these effects in the construction of our theoreti¬ 
cal predictions. 


3.2.1 Dark matter only 

The first choice we make concerns the dark matter model P^^{k), 
which is a delicate issue that has been thoroughly investigated in 
HWVH in a very similar context. Following this work, we choose 
the dark matter only model that best reproduces the results from a 
number of A^-body simulations, then implement the combined effect 
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of modified gravity, baryon feedback and massive neutrinos relative 
to this dark matter only baseline. Our dark matter only prediction is 
a hybrid model that combines the Extended Cosmic Emulator (Heit- 
mann et al. 2013) with the recalibrated HALOFIT code by Takahashi 
et al. (2012). Its convergence properties have been well examined in 
HWVH, and it was shown to have the best agreement with indepen¬ 
dent high-resolution simulation suites, compared with other models. 
In addition, HWVH examined the scatter across multiple models, and 
estimated the theoretical uncertainty on the global dark matter only 
prediction for In this paper, we also incorporate this model un¬ 
certainty in the analysis pipeline, at the level of the calculation 
(see §4.1). 


3.2.2 Neutrino and baryon feedback 

Eollowing HWVH, we model the impact of massive neutrinos and 
baryon feedback on the matter power spectrum as separate effects 
that can be expressed with multiplicative feedback terms, namely: 

p^^+’'+K-^)(^k,z) = P°^{k,z) X bli^(k,z) X hl,{k,z). (39) 


The underlying assumption is that both biases are independent, which 
is reasonable since baryons were found to have a one percent effect 
on the neutrinos bias for k < 8feMpc“^ (Bird et al. 2012). 

We compute the neutrino feedback bias term with the 
CAMB cosmological code (Lewis et al. 2000), which is reported 
to be accurate to better than 10 percent at fc = 10/iMpc“^ (Bird 
et al. 2012). We assume one massive neutrino flavour, and fix the 
cosmology at high redshift - i.e. we keep the primordial amplitude 
Aa fixed but let as vary. We justify this choice from the fact that the 
former quantity is measured very accurately by CMB observations, 
whereas our estimation of the latter quantity is much less accurate 
due to galactic and cluster bias. We construct the neutrino bias as 


r>DM+ 


{k,z) 


-Pcamb(^i 


(40) 


where the M^, (= 0.0, 0.2, 0.4, or 0.6 eV) superscript specifies the 
total neutrino mass considered, and the subscript ‘CAMB’ specifies 
that both quantities are measured from this cosmological numerical 
code. 

The baryonic feedback bias is estimated from two hydrodynam- 
ical simulations ran in the context of the Overwhelmingly Large 
(OWL) Simulation Project (Schaye et al. 2010). The dark matter only 
run (DM-ONLY) is a purely collision-less A^-body calculation and 
acts as the baseline for this baryon feedback measurement only. The 
AGN simulation run contains gas dynamics with physical prescrip¬ 
tions for cooling, heating, star formation and evolution, chemical en¬ 
richment, supernovae feedback and active galactic nuclei feedback 
(see van Daalen et al. 2011, for details about these simulations). Fol¬ 
lowing van Daalen et al. (2011); Semboloni et al. (2011), we measure 
the baryonic feedback bias by taking the ratio between the AGN and 
the DM-ONLY models^; 


P 

b'Lik,z) = 


,DM+6(m) 

OWL 


(fc, 




(41) 


where the index 6(m) refers to either DM-ONLY or AGN, and the 
subscript ‘OWL’ specifies that these quantities were measured specif¬ 
ically from the OWL simulation suite. 

Fig. 1 shows the impact of different combinations of baryons 
and massive neutrinos on the matter power spectrum. Fig. 2 and 3 



Figure 1. Combined effect from baryon feedback and massive neutrinos on 
the matter power spectrum P{k), evaluated at 2 : = 1. Results are shown with 
respect to the dark matter only non-linear predictions (thick solid line). From 
top to bottom, the (blue) dashed lines represent the effect of massive neutrinos 
with Miy = 0.2, 0.4 and O.GeV respectively. The combinations of massive 
neutrinos with baryon feedback are shown with the thin (red) solid lines. 



Figure 2. Combined effect from baryon feedback and massive neutrinos on 
the weak lensing power spectrum, assuming the source redshift distribution 
given by equation (45) and the baseline WMAP9 cosmology. As for Fig. 1, 
results of different combinations are shown with respect to the dark matter 
only non-linear predictions (thick solid line), and the sum of neutrino masses 
shown are, from top to bottom, Mu = 0.2, 0.4 and O.GeV. 




Figure 3. (left:) Combined effect from baryon feedback and massive neutri¬ 
nos on the weak lensing two-point con'elation function ^ 4 -. The open symbols 
represent our measurements from CFHTLenS data, shown with lu error bars, 
(right:) Same as the left panel, but for the estimator. We used the same y- 
axis range for both panels to emphasize on the differences across the models, 
hence the leftmost point falls outside the frame, at ^_/^5^=-3.8. 


The power spectrum measurements from the OWL simulation suite are 
publicly available at: http: //vdll. strw. leidenuniv.nl 
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show the equivalent effects on the weak lensing power spectrum Cf 
and on the shear two-point correlation function respectively. 

We can see from the three figures that all models converge to DM- 
ONLY at large scales (low k, low £ and high 0), and that the com¬ 
bined effect can suppress more than 50% of the power, depending 
on the models and neutrino mass. Also, it becomes clear that surveys 
probing small patches (restricted to i > 500 for example) would 
have difficulties to distinguish between the two feedback processes. 
This degeneracy can only be broken with the inclusion of lower I 
multipoles, where baryon feedback is minimal but massive neutrinos 
still leaves a signature (Natarajan et al. 2014). 


3.2.3 Combined feedback with modified gravity 


The evolution of perturbations in the context of large-scale structures 
has been carefully studied in f{R) and scalar-tensors theories gravity 
(Zhang 2006; Koivisto 2006; Bean et al. 2007; Song, Hu & Sawicki 
2007; Hu & Sawicki 2007; Song, Peiris, Hu 2007; Pogosian & Sil- 
vestri 2008; Carloni, Dunsby & Troisi 2008; Koyama, Taruya & Hi- 
ramatsu 2009; Motohashi, Starobinsky & Yokoyama 2009; Li & Hu 
2011; Brax et al. 2011; Li, Zhao & Koyama 2012; Linares & Mota 
2013; Taddei, Catena & Pietroni 2013; Brax & Valageas 2013). In 
this paper, we focus on the matter density power spectrum P(fc; 2 ), 
or more precisely, on the weak lensing convergence power spectrum 
Cf, which can be computed from P(fc; 2 ) through the modified Pois¬ 
son equations that relate the metric gravitational potentials to the 
matter density fluctuations. 

Therefore, before computing weak lensing statistics, we first 
need to describe gravitational clustering and the 3D matter den¬ 
sity power spectrum for all cosmological scenarios that we investi¬ 
gate. We use the approach first developed in Valageas, Nishimichi & 
Taruya (2013) for the ACDM cosmology, generalized afterwards to 
various modified-gravity scenarios in Brax & Valageas (2013, 2014). 
This is an analytical approach that combines perturbation theory up 
to one-loop order (i.e., up to order Pi, where Pl is the linear mat¬ 
ter density power spectrum) with a phenomenological halo model. 
Namely, we are splitting the matter power spectrum as: 


P{k) = P2H{k) + Piaik), 


(42) 


where P 2 H{k) is the ‘two-halo’ term associated with pairs of parti¬ 
cles that are enclosed in two different halos, whereas PiH(fc) is the 
‘one-halo’ term associated with pairs enclosed in the same halo. This 
construction allows us to obtain predictions for the non-linear mat¬ 
ter power spectrum covering the linear, quasi-linear and highly non¬ 
linear scales. We refer the reader to the work cited above for complete 
details and validations of equation (42), but nevertheless provide an 
overview of the method in the Appendix for quick reference. We note 
that other prescriptions exists for modelling P{k) in modified gravity 
scenarios, i.e. Zhao (2014) for the f{R) model. However the mod¬ 
elling we adopt here applies also to f{R) with n f 1 gravity, to 
Dilaton gravity, and in fact to any modified gravity model expressed 
in the tomographic parameterization, which makes it general and ac¬ 
curate at the same time. 

In analogy with equations (40) and (41), we define the modified 
gravity bias: 






(43) 


where MG(a) refers to the gravity model, with a= 0 corresponding 
to GR, a= [1, 2, 3, ...,15] specifying dilation models [Al, A2, A3, 
..., E4], a= [16, 17, 18] specifying f{R) models with n — 1 and 
= lO""*, 10"®, 10“®, and finally [19, 20, 21] the f{R) 
models with n = 2 and the same |/_Rg | values. The subscript ‘VNT’ 


indicates quantities that are computed in the framework of Valageas, 
Nishimichi & Taruya (2013), i.e. with equation (42). 

Bringing all the pieces together, we construct the matter power 
spectrum for any combination of baryon feedback, neutrino mass and 
modified gravity by multiplying the DM-ONLY model by the corre¬ 
sponding biases: 

7-)DM-|-iz-|-6(m)-l-MG t-jDM ;2 7,2 7 2 / aa \ 

P ^ ^ ^ = p X X X bMG(a)- (44) 

We have removed the dependences on scale and redshift for each 
of these terms to clarify the notation. This modelling assumes that 
the effect of modified gravity on the baryon and neutrino feedbacks 
can be neglected, allowing for the convenient factorization presented 
in equation (44). This seems to be a valid approximation for some 
models, as it was shown in Hammami et al. (2015) that the modified 
gravity bias measured in dark matter only matched to better than 5% 
the same measurement done in full hydrodynamical simulations, for 
/(P) models with n = 1 and I/aqI G [10“^ — 10“®]. However, 
the same group also observed larger deviations in many symmetron 
models, up to 20% by k — 10/iMpc“^ in some cases. This places a 
limit on the accuracy of equation (44), and calls for more hydrody¬ 
namical simulation runs where 6m and 6mg(c«) are merged into one 
term, 6m,MG(c«)i measured for each combination of {a, m}. This is 
unfortunately not available at the moment, hence equation (44) is cur¬ 
rently our best shot at this joint measurement. On the neutrino sector, 
results by Baldi (2014) are further encouraging: they looked at joint 
simulations of modified gravity and massive neutrinos and came to 
the conclusion that one could consider the effect of each almost in¬ 
dependently, supporting the validity of equation (44). 

For each combination, we compute predictions for the weak 
lensing quantity with equations (37) and (34). We report our results 
on P{k) and Cf in Fig. 4 and 5 respectively. Whereas modified grav¬ 
ity is generally boosting the clustering compared to a ACDM uni¬ 
verse, the inclusion of massive neutrinos and/or baryonic feedback 
is working in the opposite direction. It becomes clear that a precise 
distinction between these three feedback contributions poses a chal¬ 
lenge to clustering and weak lensing experiments. 


3.3 Data 

Our measurement of the shear correlation functions is based on 
the public release of the Canada France Hawaii Telescope Lensing 
Survey (CFHTLenS^^). The CFHTLenS covers a total area of 154 
deg^, which is reduced to 128 deg^ after masking bright stars, fore¬ 
ground moving objects and faulty CCD rows. Full details about the 
data reduction pipeline are provided in Frben et al. (2013). Source 
redshifts are obtained from the five bands u'griz photometric obser¬ 
vations (Hildebrandt et al. 2012) and were carefully tested in Ben¬ 
jamin et al. (2013); shape measurements are performed on the r- 
band images with the lensfit Bayesian code described in Miller et al. 
(2013). A detailed assessment of the residual systematics is provided 
in Heymans et al. (2012), and we refer the reader to these references 
for more information on the CFHTLenS data. 

As described in Heymans et al. (2012), the public shear data 
must be recalibrated with additive and multiplicative factors, com¬ 
monly referred to as the c and m corrections. In contrast with this ref¬ 
erence, we use a different c correction, as detailed in HWVH, which 
is less model dependent. Although the overall change on the correc¬ 
tion is marginal, the number of CFHTLenS pointings that are flagged 
as bad is reduced by almost a half. 

Following the recommendations of Heymans et al. (2012) and 


CFHTLenS: WWW. cfhtlens . org 
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Figure 4. Combined effect from baryon feedback and massive neutrinos on 
the matter power spectrum P{k) assuming different modified gravity models, 
again evaluated dX z = 1. Results are shown with respect to the dark matter 
only non-linear predictions (thick horizontal solid line). From top to bottom 
at fc = 0.2hMpc, the solid lines represent Dilaton models B4, A3, E3, D1 and 
Cl respectively. The thick red dashed lines con'espond to f{R) gravity with 
n = 1. Top to bottom are for |/hoI = 10“"^, 10“^ and 10“^ respectively. 
We do not show the n = 2 results to avoid over-crowding the figure, but 
they are qualitatively similar in shape to the n = 1 case, albeit with a smaller 
departure from ACDM. Different panels show different combinations of mas¬ 
sive neutrinos and baryon feedback on these same models, all computed with 
equation (44). 




i i 

Figure 5. Same as Fig. 4, but for the weak lensing power spectrum. 


0.14594, 0.514894, 0.498379, 0.15608, 1.74435, 0.445019, 

0.684098). There is a 0.4 percent difference in the mean redshift 
between the fit and the distribution, which yields a small error 
well below the other sources of error in our analysis. We therefore 
neglected this contribution to the systematic budget. 

We construct our shear correlation function estimator following 
Kilbinger et al. (2013): 

, , i WiWj [et{0i)et{9j) ± er{0i)er{0j)] 

iUe) = ' - ■ (46) 

WiWj 

All galaxy pairs {i,j) separated with angular distance \9i — 9j\ £9 
contribute to the same bin, with their contribution weighted by the 
product of their lensfit weights WiWj (Miller et al. 2013). The shear 
quantities et and Cr are the tangential and cross-component of the 
galaxy ellipticity, measured in the coordinate system of the galaxy 
pair. We account for the shear calibration by measuring 


1 + K{9) 


WiWj{l + mi)(l -f rrij) 
Yi,j ‘WiWj 


(47) 


and dividing by 1 + K. As a rule of thumb, AT is ~ —0.11 at all 
angular scales, with variations smaller than 0.1%. We finally exclude 
all pairs with ^ < 12 arcseconds in order to minimize contamina¬ 
tion by post stamp leakage across lensfiX. templates. We perform this 
measurement with ATHENA^®, and show our results in Fig. 3. 


3.4 Simulations 

In order to achieve a high precision cosmic shear measurement, not 
only must the data be thoroughly tested for subtle systematics resid¬ 
uals, but the sampling variance must be accurately estimated, a quan¬ 
tity that is very hard to assess from the data. To overcome this dif¬ 
ficulty, we rely on a suite of weak lensing simulations based on 
WMAP9 -I- SN -I- BAO cosmology. As detailed in Harnois-Deraps & 
van Waerbeke (2015), the SLICS-LE suite consists of 60 sq. degrees 
light cones extracted from 500 independent A^-body realizations. The 
numerical weak lensing signal is precise to better than 10 percent for 
with 9 > 0.4' (and 9 > 5' for ^_). We construct the mock maps 
by combining the different redshift planes with a redshift source dis¬ 
tribution that mimics that of the data. We then sample the simulated 
shear maps with 10® points randomly located, and compute the shear 
two-point correlation functions of these mock ‘galaxies’ with the 
same pipeline as the data (i.e. from equation 46). 


4 RESULTS 

In this Section, we first review our error budget, we then describe how 
different components combine in our model rejection procedure, and 
finally present our results. 


Benjamin et al. (2013), we minimize the systematic contamination 
from badly reconstructed photometric redshifts by applying the se¬ 
lection cut 0.4 < Zphot < 1.3. We construct the redshift distribu¬ 
tion n{z) for the selected galaxies from the Zeni-fit-weighted stacked 
probability distribution functions of the galaxy sample. As shown in 
HWMH, the distribution is well described by the following analytical 
expression: 

niz) = -f 

where (No, zo,cro, Ni, zi,ai, N 2 , Z 2 ,cr 2 ) = (0.14438, 0.760574, 


4.1 Error budget 

This analysis closely follows that of HWVH; we summarize here the 
main ingredients, and refer the reader to the reference for more de¬ 
tails. The sources of error in this analysis can be broken into three 
terms: 1- uncertainty on the cosmic shear measurement, 2- uncer¬ 
tainty in the theoretical model describing the non-linear regime of 
structure formation, and 3- uncertainty on the fiducial cosmology. 

1- The error on our cosmic shear measurement is dominated 
by shape noise at small angles and sampling variance at large angles. 

ATHENA: http : / /cosmostat. org/athena . html 
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Figure 6. Same as Fig. 4, but for 5+- The open symbols represent our mea¬ 
surements from the CFHTLenS data, exactly as in the right panel of Fig. 3. 
Shown are the Icr error bars. 




^[arcmin] ^[arcmin] 

Figure 7. Same as Fig. 6, but for . Note the different ty-axis range compared 
to Fig. 6. 

The angular scales at which these two errors contribute equally occur 
at 0 — 2 and 30 arc minutes for and respectively. In addition, 
the variance-shape noise mixed term contributes to about a third of 
the error on at large angles, but is negligible in , as seen in Kil- 
binger et al. (2013). We have estimated the sampling variance from 
the SLICS-LE weak lensing simulations, and added an extra contri¬ 
bution from the halo sampling variance, following the modelling of 
Sato et al. (2009), which provides at most a 10 percent correction 
on the overall error. Our measurement is minimally affected by in¬ 
trinsic alignment of galaxies, since we do not perform tomographic 
analysis (see Heymans et al. 2013, for more details on intrinsic align¬ 
ments in the CFHTLenS data). The error from shape reconstruction 
is already included in the statistical term, hence does not require an 
extra term. Photometric redshift uncertainty enters the measurement 
through modification of the source distribution n(z), but this effect 
is negligible compared with other sources of error hence is not in¬ 
cluded. 

2- The uncertainty on the dark matter only non-linear model has 
been carefully assessed in HWVH by comparing five different pre¬ 
dictions: F1ALOFIT2012, F1ALOFIT2011 -I- small scale empirical re¬ 


calibration, Cosmic Emulator -l- power law graft. Cosmic Emulator -l- 
F1ALOFIT2012 graft and, finally, the mean over five independent high 
resolution simulations - the SLICS-HR suite described in Hamois- 
Deraps & van Waerbeke (2015). These five models agree very well 
over most angular scales, and the Icr scatter among them is taken as 
the (^-dependent) theoretical error. 

3- The uncertainty in the cosmological parameters is set by the 
WMAP9 precision (Hinshaw et al. 2013), whose dominant contribu¬ 
tion on the weak lensing uncertainty arises via the parameters Hm 
and As- With the inclusion of the BAO and SN external data, these 
two parameters are allowed a 3.4 and 3.3 percent variation about 
their mean values (Icr). Since the amplitude of roughly scales as 
we expect the combined error to be of the order 5 percent 
of the ACDM baseline signal, assuming no prior on the joint contour. 

Note that the cosmological error and the modelling error (terms 
2 and 3) enter in our analysis as systematic uncertainties, therefore 
we add them in quadrature and marginalize over them (see details in 
§4.2). Also note that the Planck cosmology {Hm, As} falls within 
our 3cr search limits, although closer to the edge of the search zone. 

4.2 Model rejection strategy 

As seen in Fig. 6 and 7, the effects of baryons, massive neutrinos 
and modified gravity are significantly degenerate on Given the 
noise levels in the current data and the number of internal parameters 
that describe these different mechanisms, performing a full MCMC 
analysis is not convenient to extract meaningful constraints. A more 
appropriate and direct way is to sample a finite set of model combi¬ 
nations and examine their agreement with the data. This case-by-case 
strategy has the potential to reject models that are inconsistent with 
the data, which can then be translated into constraints on the under¬ 
lying free parameters. 

The metric we adopt for this type of analysis is the p-value, 
which measures the probability that the data is consistent with the 
model, if the model is true. It is given by the integral over the 
probability density function, where the lower bound is the measured 
X^ and the upper bound is infinity. As a rule of thumb, models with 
p-values < 10% are rejected with more than 90% confidence, and 
Icr, 2 (t, 3(7... rejection measurements are obtained for p-values of 
0.317, 0.046, 0.003.... Our strategy therefore consists to measure the 
X^ and p-value associated with each combination of baryon feed¬ 
back, neutrino mass and modified gravity model, and to flag every 
combination withp < 0.1 as being disfavoured. 

The uncertainty arising from statistical and sampling variance 
naturally enters this calculation through the evaluation of the x^. 
which involves the inversion of the cosmic shear covariance matrix. 
The systematic uncertainty, however, is trickier to capture. In our cos¬ 
mic shear measurement, it mainly manifests itself as shifts in the am¬ 
plitude of the signal, as described in §4.1. The systematic error is 
higher at smaller angles and represents at most an error of ~ 9% on 
the model amplitude. In order to marginalize over this effect, for 
each model, we allow the theoretical signal to shift up and down by 
3crsyst, corresponding to vertical excursion of 27% in Fig. 6 and 7, 
keeping the error bars (statistical -l- sampling) fixed. We then com¬ 
pute an array of p-values in this excursion range, and record only the 
largest measurement (i.e. the least restrictive). 

The exact number of degrees of freedom (d.o.f.) that enters the 
X^ distribution function must be carefully chosen. To begin with, 
each of the two cosmic shear observables is organized in 11 angular 
bins, yielding a maximum of 22 d.o.f. However, assigning one d.o.f. 
per angular bin would be incorrect, for the following reason. In a sta¬ 
tistical sense, our model rejection method is completely equivalent 
to fitting the parameter combination (AjH}/) from the amplitude of 
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the signals, followed by an extraction of the most likely neutrino 
mass for each baryon feedback and modified gravity model from the 
largest p-value. This implies that the number of degrees of freedom 
should be reduced by two (one for fitting one for fitting M„) 

in the conversion between and p-values. 

Note that for a given angular scale, both and probe differ¬ 
ent physical scales, the latter focusing on structures about five times 
smaller. It is therefore relevant to examine the constraining power 
of first, and to add to the data vector as a second step. When 
both are combined, the full data covariance matrix involves the cross¬ 
correlation region, as described in HWVH. The resulting p-values are 
summarized for all our results in Table 2, for ^ 0.2 eV. No con¬ 
clusions can be drawn from models with higher total neutrino masses, 
as the p-values for any combination is always greater than 0.175. The 
models rejected at more than l.fiTcr (i.e. 90% Cl) are highlighted in 
bold font. 

4.3 Discussion 

One of the main results recovered from Table 2 is that the f{R) 
model with {|/ho|,m} = {10““^, 1} is strongly disfavoured by the 
cosmic shear data, regardless of the baryonic feedback model or sum 
of neutrino mass, consistent with independent constraints. The f{R) 
and f{R) + AGN models are rejected by at least 3a, but combina¬ 
tions including massive (0.2eV) neutrinos tend to weaken this con¬ 
straints. This can be understood by the fact that massive neutrinos 
and modified gravity partly compensate for one another, reducing the 
global departure from ACDM. Also, f{R, n=2) models are generally 
in better agreement with the data compared to their (n=l) counter¬ 
part. This is so simply because higher values of n rapidly suppress 
the f(R) term, hence deviations from GR, as seen in equation (22). 

The next important result is that the rejection of massless neu¬ 
trinos -I- DM-ONLY is robust against all modified gravity models 
we have tested, and typically made stronger. The cosmic shear data 
clearly prefers lower values of at small angular scales, and modi¬ 
fied gravity pulls the other way. 

The inclusion of baryon feedback reduces to about two thirds 
the number of models rejected with 90% Cl. For instance, dilation 
models A2, A3, B3, B4, Cl, C3, Dl, E3 and E4 are disfavoured; 
these are the most discrepant with GR-tACDM. Referring to Table 1 
and the model descriptions in §2.1, this can be interpreted as follow. 
In a tomographic parameterization of modified gravity centred on 
{mo, r,/3o, s} = {0.334, 1.0, 0.5, 0.24}, excursion in the s, /3o and 
r directions are studied with models A, B and C respectively, and 
the data favours lower parameters values. Model E explores the mg 
direction, and the data prefers higher values. Model D explores the 
diagonal direction in the {mo, s} plane at fixed A 2 (see equation 8), 
where here we observe instead that the data prefers lower mo values. 

We note that there is a mild effect seen in the ‘AGN’ column of 
Table 2, where the addition of to the data vector sometimes in¬ 
creases the p-value by a small amount. This can be attributed to the 
fact that at small angles, prefers amplitude even lower than ^+, 
compared to the DM-ONLY model. Adding baryon feedback there¬ 
fore produces a lower rejection rate in the former than in the latter 
quantity. 

When neutrino masses are allowed to reach 0.2eV, only the 
f(R){10~‘^, 1}, 2} and the Dilaton B4 and E4 models 

remain in tension with the data. With AGN - 1 - Mv=Q-2eN, no mod¬ 
els are rejected, aside from the most extreme case considered in this 
paper: /(i?){10"'‘, 1}. 

This means that given the current cosmic shear data and levels 
of systematics, it is possible to accommodate most models, as long 
as either massive neutrinos or strong baryon feedback mechanisms 


counter-balance the effect of the fifth force on the matter clustering. 
As upcoming independent cosmological probes will tighten the un¬ 
certainty on neutrino masses and significantly improve the statistical 
and sampling errors, we expect the next generation of such analysis 
to be much more constraining. Once at this stage, it will be instruc¬ 
tive to propagate our measurements onto {m(a), /3(a)} contours and 
provide a Fisher matrix for joint probes analyses. If, for instance, the 
total mass turns out to be much smaller than 0.2eV, then the current 
AGN column should gives a very good approximation of the rejection 
power from the CFHTLenS cosmic shear data. Precise modelling of 
the baryon feedback is likely to take more time to reach, due to the 
higher level of complexity intrinsic to these astrophysical phenom¬ 
ena. Intermediate solutions will involve a series of tuneable parame¬ 
ters, also to be constrained. 

On a separate note, we stress that the constraints can be fur¬ 
ther tightened using additional information about the weak-lensing 
observables, such as the non-Gaussian features (Munshi, van Waer- 
beke, Smidt & Coles 2012), or by combining the results with external 
probes such as redshift distortions, peculiar velocity, etc. 


5 CONCLUSION 

Cosmic shear is a promising tool for probing deviations from GR, 
since these are maximal at scales of a few Mpc, where the lensing 
signal-to-noise ratio is the highest. These same scales are very chal¬ 
lenging to probe with other types of large scale structure observables, 
mainly because of the galaxy bias that is largely unknown. At the 
same time, this complimentarity offers a number of opportunities for 
strong constraints based on joint data sets. 

One of the main challenge in working with these non-linear 
scales is the large theoretical uncertainties due to the unknown neu¬ 
trino masses, the precise baryonic feedback mechanisms and, to a 
lesser extend, inaccuracies in the clustering of dark matter. However, 
a lot of effort is invested in all these areas, such that it becomes pos¬ 
sible to place joint constraints on these degenerate physical effects. 

This paper presents the first constraints on modified gravity ob¬ 
tained from cosmic shear measurements alone; the results are derived 
by studying the impact of modified gravity on matter clustering and 
comparing the predictions with the public CFHTLenS data. Limiting 
the background ACDM cosmology to the 30- range in {As, Dm} al¬ 
lowed by WMAP9 + SN + BAG, we compared the data against 
predictions including f{R) and Dilaton models in a number of pa¬ 
rameter configurations. We carried a model rejection analysis ac¬ 
counting for possible degeneracies with massive neutrinos and bary¬ 
onic feedback mechanisms, and investigated which combinations of 
models were mostly disfavoured by the data. As summarized in Table 
2, the f{R) model with I/hqI = 10“"* is strongly disfavoured even 
in the presence of realistic levels of baryonic feedback and massive 
neutrinos reaching M„=0.2eV. A universe with no baryonic feed¬ 
back and massless neutrinos is also rejected with 2a or above in 
most modified gravity scenarios. We are not yet able to identify a 
preferred model with the current level of statistical accuracy, but we 
expect future weak lensing experiments to improve significantly in 
this direction. 

In our analyses, we have use the simplifying assumption that 
the biases due to massive neutrinos, baryon feedback and modified 
gravity were uncorrelated, which is justified to some extend based 
on the several numerical results. However, precise correlations will 
need to be studied for a number of models, a task that involves suites 
of large cosmological hydrodynamical simulations including all the 
ingredients at once. 

One important future tasks will be to map observational detec- 
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Table 2. Distribution of p-values for different combination of baryon feedback models, neutrino masses and gravity models (see main text for details). The 
parameters listed in the leftmost column of the f{R) models are {n, I/aq |}. Dilaton models are described in Table 1. For this calculation, we fit all data in the 
range 0.2 < 9 < 167 arcminutes. Specifically, each entry in this Table represents the largest p-value probed inside a Susyst region about the mean of the model. 
Values in bold face highlight the model combinations that are excluded by the data with more than 1.64 (t significance (p-value < 0.1, equivalent to a confidence 
interval (Cl) of 90%). Models with > 0.2 eV are not listed, as none has value lower than 0.176. 


Model 

DM-ONLY 

AGN 

0.2eV 

«+ S+C- 

AGN+0.2eV 

General Relativity 

ACDM 

0.132 

0.065 

0.119 

0.150 

0.331 

0.297 

0.289 

0.444 

Generalized Dilaton 

Al 

0.126 

0.058 

0.116 

0.141 

0.323 

0.282 

0.284 

0.431 

A2 

0.088 

0.030 

0.093 

0.099 

0.269 

0.203 

0.256 

0.363 

A3 

0.037 

0.008 

0.060 

0.051 

0.184 

0.105 

0.207 

0.265 

B1 

0.120 

0.054 

0.113 

0.135 

0.315 

0.273 

0.281 

0.424 

B3 

0.043 

0.010 

0.064 

0.054 

0.195 

0.113 

0.212 

0.272 

B4 

0.008 

0.001 

0.032 

0.019 

0.107 

0.037 

0.150 

0.161 

Cl 

0.080 

0.022 

0.100 

0.098 

0.259 

0.171 

0.271 

0.365 

C3 

0.104 

0.040 

0.098 

0.111 

0.293 

0.239 

0.261 

0.385 

C4 

0.114 

0.049 

0.104 

0.122 

0.307 

0.259 

0.267 

0.402 

D1 

0.063 

0.013 

0.100 

0.085 

0.232 

0.132 

0.276 

0.347 

D3 

0.127 

0.060 

0.115 

0.141 

0.325 

0.286 

0.282 

0.431 

D4 

0.131 

0.064 

0.119 

0.149 

0.331 

0.295 

0.288 

0.441 

El 

0.118 

0.049 

0.117 

0.135 

0.314 

0.260 

0.289 

0.425 

E3 

0.047 

0.013 

0.053 

0.049 

0.200 

0.132 

0.186 

0.257 

E4 

0.026 

0.007 

0.032 

0.027 

0.156 

0.094 

0.138 

0.188 

f{R) 

{1,10-^} 

0.001 

0.000 

0.005 

0.003 

0.051 

0.018 

0.054 

0.057 

{1,10-5} 

0.058 

0.013 

0.072 

0.062 

0.222 

0.134 

0.222 

0.292 

{1,10-®} 

0.129 

0.054 

0.125 

0.145 

0.328 

0.271 

0.298 

0.437 

to 

o 

1 

0.011 

0.003 

0.020 

0.014 

0.112 

0.056 

0.104 

0.131 

(2,10-5} 

0.095 

0.030 

0.094 

0.097 

0.277 

0.200 

0.254 

0.354 

{2,10-®} 

0.137 

0.063 

0.126 

0.154 

0.338 

0.292 

0.299 

0.449 


tions of modifications to GR onto parameter constraints such as the 
{m{a), P{a)} pair. However, the current data is not quite there yet. 
Several theories can accommodate similar phenomenological effects, 
and model-independent parameterizations such as that presented in 
Leonard, Baker & Ferreira (2015) might prove helpful for this. 

This paper used the impact of modified gravity on the clustering 
properties of matter and their propagation onto the weak lensing cos¬ 
mic shear signal. Other avenues of probing deviations from GR with 
weak lensing data are complimentary, including direct combinations 
with baryonic probes or tomographic decomposition, and worth ex¬ 
ploring in the near future. 
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APPENDIX: Details on the Theoretical Modelling of P(fc) in 
Modified Gravity Scenarios 

This Appendix discusses the construction strategy of the matter den¬ 
sity power spectrum P(fc; 2 ) in the presence of f{R) or Dilaton mod¬ 
ifications to General Relativity; full details are provided in the refer¬ 
ences contained herein. 

1 Two-halo term: P 2 H(fc) 

The power spectrum analytical prediction is constructed from a halo 
model approach, following equation (42). The two-halo term, which 
dominates on large scales, is computed from a Lagrangian-space re¬ 
summation of standard perturbation theory that is exact up to order 
Pl and contains partial resummations of higher order terms. It is 
also supplemented with non-perturbative contributions that take into 


© 0000 RAS, MNRAS 000, 000-000 


Testing Modified Gravity with Cosmic Shear 13 


account some aspects of shell crossing and ensure that all particle 
pairs are counted only once in the sum. Within this framework, the 
large-scale term P 2 H if) essentially contains no free parameter. It can 
therefore be computed in ACDM and modified-gravity scenarios by 
using perturbation theory up to order Pi (which requires going to 
order in terms of the density field itself). 

In the case of the ACDM cosmology, this perturbative ex¬ 
pansion follows the standard approach (Bernardeau, Colombi, Gaz- 
tanaga & Scoccimarro 2002), where the density and velocity fields 
are written as perturbative expansions over powers of the linear den¬ 
sity field Sl', subsequent orders are computed by substituting into the 
continuity and Euler-Poisson equations. In the case of the modified- 
gravity scenarios considered in this paper, we follow the same ap¬ 
proach but require an additional expansion to write the fifth force in 
terms of the non-linear density fluctuations. Indeed, using the quasi¬ 
static approximation, we can relate the scalar field ip to the matter 
density field p, typically through a non-linear Klein-Gordon equa¬ 
tion. Then, we can solve for p as an expansion over the non-linear 
density fluctuations Sp. This allows us to obtain both the Newtonian 
potential and the fifth-force potential as functionals of the non-linear 
matter density fluctuations. However, while the Newtonian potential 
is given by the linear Poisson equation, the fifth-force potential is 
usually given by a non-linear equation that involves new time and 
scale dependences. In terms of the diagrammatic expansion of the 
non-linear power spectrum P{k) over PL{k), this implies that the 
linear propagators and the vertices are modified with new diagrams 
associated with the new non-linearity of the modified Poisson equa¬ 
tion. See Brax & Valageas (2013) for more explanations. 

2 One-halo term: PiH(fc) 

The one-halo term is obtained from the halo mass function and the 
halo density profile, with the addition of a counter-term first intro¬ 
duced in Valageas & Nishimichi (2011) that arises from mass conser¬ 
vation. This also ensures that Pmik) decays at low k and becomes 
subdominant as compared with P 2 ii{k), whereas the usual formula¬ 
tion gives a spurious white-noise tail that dominates on very large 
scales. We take into account the impact of modified gravity through 
its effect on the halo mass function (i.e., through the acceleration or 
slowing down of the spherical collapse), but neglect the impact on the 
halo shape and profile. This should be sufficient for our purposes, be¬ 
cause we only consider cosmologies that remain close to the ACDM 
reference, and these modified gravity models have a much stronger 
impact on the halo mass function, especially on its large-mass tail, 
than on the halo profile. As shown in Valageas (2013), at z = 0, 
a 10% change to the mass-concentration relation only yields a 2% 
change of P{k) at l/iMpc~^, whereas a 10% change to the halo 
mass function yields a 2% change of P{k) at 0.35h,Mpc~^ and a 
7.5% change at l/iMpc”^. Generally, the concentration parameter 
always remains in the range 3—10 for typical halos and does not 
vary by much more than 10% for realistic scenarios, whereas the 
mass function at M ~ 5 x Mq can vary by more than 50% 

(Lombriser, Koyama, Zhao & Li 2012; Lombriser, Li, Koyama & 
Zhao 2013). The interior of haloes are mostly affected by screening 
anyway, further justifying this approximation. 

3 Comparison with numerical simulations 

The modelling described above for the matter density power spec¬ 
trum has been checked in details against numerical simulations in 
Valageas, Nishimichi & Taruya (2013) for ACDM cosmologies, and 
in Brax & Valageas (2013) for the class of modified gravity mod¬ 
els that we consider in this paper. In the case of ACDM, it provides 


an accuracy of 2% up to comoving wavenumber k ~ 0.9/iMpc“^, 
and 15% up to fc = 15/iMpc“^ down to z = 0.35. In terms of the 
real-space correlation function, this translates into an accuracy of 5% 
down to the comoving scale r = 0.15/i“^Mpc. For the f{R) theo¬ 
ries, it is able to reproduce very well the deviations from the ACDM 
scenarios up to fc = 3h,Mpc“^ (the highest wave number available 
from the simulations) at z = 0, for | = lO”'^, 10“®, and 10“®. 
In particular, it accurately captures the relative effect on the power 
compared to the ACDM reference due to the non-linear Chameleon 
mechanism. For the Dilaton models, the agreement with the numer¬ 
ical simulations depends somewhat on the model parameters but it 
typically gives a good quantitative estimate of the deviations from 
ACDM up to fc = 5/iMpc“^ (the highest wave number available 
from the simulations). When there is a noticeable departure from the 
simulations, it corresponds to an underestimation of the amplification 
of the power spectrum at fc > 2/iMpc“^, which may be due to our 
neglect of the impact of modified gravity on the halo concentration 
parameter. Therefore, in such cases our approach provides a conser¬ 
vative estimate of the deviations from ACDM. Again, this modelling 
is able to capture the decrease of the deviations from the ACDM ref¬ 
erence due to the non-linear Damour-Polyakov mechanism. 
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